Origin of second-harmonic generation enhancement in optical split-ring resonators 
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We present a study of the second-order nonlinear optical properties of metal-based metamaterials. 
A hydrodynamic model for electronic response is used, in which nonlinear surface contributions are 
expressed in terms of the bulk polarization. The model is in good agreement with published exper- 
imental results, and clarifies the mechanisms contributing to the nonlinear response. In particular, 
we show that the reported enhancement of second-harmonic in split-ring resonator based media is 
driven by the electric rather than the magnetic properties of the structure. 
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Metamaterials (MMs) are artificially structured media 
whose collective electromagnetic properties derive from 
the geometry of sub-wavelength inclusions. To date, the 
most common MM designs have made use of inclusions 
formed by conducting materials that function as sub- 
wavelength electrical circuits. These conductor based 
MMs have proven adept at mimicking a wide variety of 
linear electromagnetic responses, providing a new venue 
to explore otherwise inaccessible concepts^. In the con- 
text of nonlinear response, however, artificial materials 
may offer even greater opportunities, due to the inher- 
ently inhomogeneous local field distribution that exists 
within and around MM inclusions. By carefully struc- 
turing the inclusion geometry, extremely large field en- 
hancement regions can be produced that can dominate 
and enhance the effective nonlinear response of the com- 
posite. 

The enhancement of nonlinear processes by MMs has 
been demonstrated at radio and microwave frequencies, 
using packaged components, such as varactor diodes, 
to introduce nonlinearity into the gaps of metal MM 
inclusionJ^. However, to achieve nonlinear optical ma- 
terials at higher wavelengths, a simple scaling of these 
prototype structures to higher frequencies (e.g., beyond 
a few terahertz) will not suffice. First, the response of 
most metals changes from conductor-like to dielectric-like 
at frequencies above a few terahertz, with absorption in- 
creasing significantly as the fields are able to penetrate 
further into the metal. Second, packaged semiconductor 
components are not readily available at frequencies above 
100 GHz. 

While metals and conductors may possess undesirable 
properties at optical wavelengths, such as increased ab- 
sorption, they also possess unique and potentially advan- 
tageous properties. In addition to large field enhance- 
ments, metal nanostructures also support intrinsic non- 
linearities that relate to the dynamics of free and bound 
charge carriers. As a result, metals possess some of the 
largest nonlinear susceptibilities known. Examples in- 
clude the large x*^^^ values of gold or silver, for example, 
suggesting that metals can serve both to form the linear 
MM response by tailoring the structure, while serving as 
the source of nonlinearity for nonlinear optical MMs. 



The second-order nonlinearity in metals arises from 
both volume and surface contributions. Nonlinear sur- 
face contributions are strictly related to the response of 
the electrons within the Thomas-Fermi screening length 
{\tf ^0.1 nm for gold) from the metal boundaries. In 
this sub-nanometer realm, electron-electron interactions 
become appreciable and non-local effects must be taken 
into account. Moreover, since metals are centrosymmet- 
ric, they do not possess an inherent x^^'' nonlinearity. 
However, the surface of a metal can break spatial sym- 
metry and provide a mechanism for an effective x*-^-* non- 
linearity. This homogenized x^"^^ nonlinear response be- 
comes highly dependent on the metal geometry, making 
it inherently a MM construct. 

A steady stream of works concerning harmonic gener- 
ation from metallic nano-structures has been published 
recentljJsHZl Second-harmonic generation (SHG) has also 
been observed experimentally from a variety of metal 
nanoparticle system^^HHI^ and specifically from optical 
split-ring resonators (SRRs), which at first app eared to 
exhibit an anomalously high conversion efficiencj^i^IIIIlol 
The magnetic resonances associate with SRRs have 
raised speculation that the enhanced SHG results from a 
strong nonlinear magnetic response associated with the 
Lorenz forc^^^. Yet, a convincing explanation of the na- 
ture of this conversion enhancement remains lacking in 
the literature. 

In this Letter, we show that the basic characteristics 
of nonlinear optical SRRs may be explained solely by the 
electric properties of the structure rather than its mag- 
netic response, in co ntrast to the conclusions drawn from 
previous work^^^li^l xhe nonlinear optical response of 
the charge carriers in the metal here is described by a hy- 
drodynamic model, which includes the effects of pressure 
associated with the electron gas. To facilitate the nu- 
merical models and remove ambiguities associated with 
additional boundary conditions, we develop an expres- 
sion for the effective second-harmonic surface currents in 
terms of the polarization vector in the bulk regions. This 
expression allows us to easily study the SHG process from 
metal particles of arbitrary shape. 

In the context of the hydrodynamic model, the elec- 
tron fiuid density, n(r,t), and the electron velocity field. 
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v(r,t), satisfy Euler's equation: 

— + (vV)v + 7v=— (E + vxB)-^Vn, (1) 
ot m* n 

along with the continuity equation, V • J = —en, with 
J = env. In Eq. 0, m* is the effective electron mass, 
and 7 is the electron collision rate. The last term in 
the equation is due to the electron gas pressure, here 
described within the Thomas-Fermi model, with /3 pro- 
portional to the Fermi velocity vp ■ Combining the con- 
tinuity equation with Eq. ([T]), and expanding all fields 
in a perturbative manner one finds that the free electron 
polarization P = J satisfies the following set of inhomo- 
geneous equation^! 



/3'V(V-Pi) 



zc.7)Pi = -^Ei, (2a) 



/J^V (V • P2) + (4lj2 + 2iw7) P2 = -^E2 + Snl, 



where the nonlinear source, Snl, is given by: 



-^[(V-Pi)Pi 



r?Pi X Bi 
KPi-V)Pi] 



(2b) 



(3) 



The subscripts refer to fundamental and second- 
harmonic fields respectively. In deriving Eqs. ^ a har- 
monic time dependence (e~''^*) has been assumed. To- 
gether with Maxwell's equations, these equations de- 
scribe the fundamental and second-harmonic polarization 
vectors, respectively, and hold under the assumption that 
the fundamental field remains undepleted. The nonlinear 
source given by Eq. ([3| groups surface and bulk second- 
harmonic contributions. Specifically we have the non- 
linear Coulomb term, (referred to as a quadrupole-like 
term by virtue of its form) proportional to Ei (V • Pi), 
the magnetic Lorentz force contribution. Pi x Bi, and 
convective terms (V • Pi)Pi and (Pi • V)Pi. 

The effect of the electron gas pressure on the polar- 
ization vectors Pi and P2 in Eqs. ([2]) is a linear, nonlo- 
cal contribution of the form /3^V (V • P). This term has 
been predicted to be responsible for unus ual, r esonant - 
like phenomena above the plasma frequencjl^lEll. in gen- 
eral, this term becomes important in a region of or- 
der of Xtf in the vicinity of the metal surface, where 
electron-electron interactions dominate the nonlinear 
surface sources given by Eq. ([S]). However, the linear 



non-local term may be safely neglected in Eq. ( 2b ) , since 
it does not affect the amount of generated harmonicP^. 

The presence of spatial derivatives (nonlocality) in 
the description of the fundamental polarization vector of 
Eq. (2a) requires the specifications of additional bound- 
ary condition s to solve the electromagnetic boundary 
value problerrPSmi, From the continuity equation and 
Gauss' theorem one obtains that n • P = at the 
boundary^^ A consequence of the linear pressure term 



is the non-zero extension of the induced surface charge 
inside the metal. That is, the induced electron charge 
density, • Pi , is zero in the bulk region, and rapidly 
changes near the metal surface, where it reaches its max- 
imum value. This behavior reflects the variation expe- 
rienced by the normal component of Pi, which continu- 
ously goes to zero at the metal boundary. 

The hydrodynamic model is a simplistic model of elec- 
trons that nevertheless gives a fairly accurate description 
of linear and nonlinear processes occurring at the sur- 
face of metallic structures. However, the simultaneous 
manifestation of microscopic and macroscopic scales in 
the problem makes the resolution of the system of (non- 
local and nonlinear) equations quite complex, and ordi- 
narily necessitates considerable computational resources 
even for particles whose dimensions are a few tens of 
nanometers. Deep inside the metal, in the bulk region, 
the electron pressure may be neglected, since it leads 
to corrections of order (At_f/A)^ <C l^'^. On the other 
hand, second-harmonic conversion efficiency strongly de- 
pends on the behavior of Pi at the metal surface where 
the impact of the pressure term may be critical. More- 
over, if the electron pressure in Eq. ( 2a I is neglected at 
the surface, the theory becomes inherently ambiguous, as 
pointed out by Sipe et alP^. 

To understand where the ambiguity arises consider the 
nonlinear Coulomb term in Eq. ([s]) , which is proportional 
to (V • Pi)Ei. In the free electron limit (/3 — 0) the nor- 
mal component of the electric field En = n • Ei changes 
discontinuously across the metal boundaries. This dis- 
continuity is attributed to charge accumulation occurring 
at the surface, which assumes the form of a Dirac delta 
function, namely, V • Pi = P„(5o(n • r), where P„ is the 
the bulk polarization component normal to the surface. 
It follows that the quantity (V • Pi)£'„ ~ PnSoEn is not 
well defined and cannot be integrated. The same analy- 
sis applies to the convective terms. In early works this 
problem was circumvented by introducing phenomeno- 
logical coefficients of order unitj-^^, or through an effec- 
tive plasma frequencj'^S^ that incorporates the details of 
the charge distribution near the surface, an effect that 
was neglected in Ref. However, an analysis of the ef- 
fect of pressure on the amount of converted energy shows 
that the total SHG converges to an asymptotic value as 
/3 tends toward zero. We verified that for typical values 
of /3 (of order of 10^ m/s ), the actual converted energy 
differs by not more than 3% from its asymptotic value. 
Exploring the limit for /3 — > seems then a very good way 
to get an approximate solution without having to solve 
the complex nonlocal equations, as discussed below. 

In this limit the surface layer is so small compared to 
the size of the nanoparticle that only the derivative along 
the direction, ^ = r • n, normal to the metal boundary 
matters. Let us choose the metal boundary such that for 
^ < the electron pressure may be neglected (/3 = 0); the 
region < ^ < I, where I ~ Xtf, represents the surface 
layer where the pressure is important. In this region, the 
parallel component of the vector Pi is nearly constant. 
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p|'(^) ~ -F'i(0 ), while its normal component may 
be written as: 



(4) 



where (t(^) is an unknown, rapidly varying function that 
goes from cr(0) = 1 to cr(/) = (as imposed by boundary 
conditions). It is straightforward to show that starting 
with Eq. (2b), that the nonlinear surface polarization is 
given by: 
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(Pi •V)Pi 



(5) 



where the second-harmonic pressure term, /3'^V(V • P2), 
has been neglected, as already pointed out. In writing 
Eq. ([5]) we have used Eq. (2a) to express the electric 
field El as a function of Pi. The crucial quantity here 



is the effective nonlinear surface current density, defined 
as: K^TL = 2iuj J^^Jq' Pf (r)dr, where the integral is per- 
formed across the surface layer. Using Eq. ([5| we obtain: 
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where /q tr'd^ = — 1 and J^^aa'd^ = —1/2 have been 
used, with the prime representing the derivative with re- 
spect to ^, and the unit vector t pointing in the direction 
n X Pf . Eq. ^ states that surface contributions to the 
second-harmonic polarization may be approximated by 
an effective nonlinear current sheet at the surface of the 
nanoparticle. These currents are related to the polar- 
ization values in the bulk region and do not require the 
resolution of nonlocal equations. That is, the volume 
sources are calculated by assuming /3 — while the sur- 
face nonlinear currents are given by our Eq. ^ . In the 
limit of validity, this approach provides a description of 
the SHG process that may be easily implemented in full 
3D simulations. We tested this approach in the case of 
SHG from a silver slab and compared our results with 
experimental dataP^. The results obtained using Eq. ^ 
agreed very well with both the solution given by the com- 
plete Eqs. (2| and the experimental data. 

3D metal nanoparticle systems are of considerable in- 
terest for nonlinear media, because both the field en- 
hancement regions and the surface morphology strongly 
impact SHG. In one of the first experiments of its kind, 
planar arrays of varying MM inclusions were shown 
to produce second-harmonic light, with efficiencies that 
varied according to the shape of the nanostructured 
inclusions'^. To investigate the mechanism of SHG 
from metal nanoparticles, we choose the same types of 
nanoparticles studied in Ref. [12^. Our results are sum- 
marized in Fig. [T] for the different kinds of nanoparti- 
cles. To ease the computational burden and to simplify 
the system studied, we assumed the nanoparticles to be 
surrounded by air rather than including the substrate 
material used in the experiments. We do not expect sig- 
nificant changes in the observed mechanism of SHG due 
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FIG. 1. (Color on line) Second-harmonic conversion efficiency 
for different 3D gold nanoparticles. All the particles are 
20 nm thick. The linear transmittance (solid line) and re- 
flectance (dashed line) at normal incidence are shown in the 
second column for the variety of nanoparticles. The verti- 
cal lines indicate the fundamental (dotted red line), and the 
second-harmonic (blue dashed-dot line) wavelength, respec- 
tively. In the third column, the arrows denote the polarization 
state. The fundamental field is represented in red (light gray) , 
the second-harmonic field is blue (gray). The relative con- 
version efficiencies normalized with respect to the U-shaped 
nanoparticle efficiency are also reported (the corresponding 
experiment of Ref.^are shown inside brackets). The follow- 
ing values for the parameters have been used: m* = me, 
no = 5.7 X 10^^ cm-^ and 7 = 1.07 x lO" s"\ 



to the exclusion of the substrate material. The structure 
studied extends periodically in the x and y directions, 
so that only a single unit cell is needed in the computa- 
tional space. To avoid possible numerical artifacts due 
to the field localization near metal corners, we consid- 
ered rounded corner geometries with a radius of curva- 
ture of 5 nm. The geometrical parameters were chosen so 
that the nanoparticles would display a resonance around 
ApF = 1.5 /im, to which the fundamental field is tuned. 
The conversion efficiencies, 77, assume an average pump 
intensity of ~ 55 MW/cm^ with the electric pumping 
field polarized along the x-direction. 

We find qualitatively good agreement with the ex- 
perimental data of Ref. [ [12] , where second- and third- 
harmonic generation were experimentally investigated for 
rectangular, T- and U-shaped gold nanoparticles respec- 
tively. We find the relative efficiencies for normal inci- 
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dence normalized with respect to the U-shaped nanopar- 
ticle efBciency to be about 1.7% for T-shaped nanopar- 
ticles (compared with 2.1% for the experimental data), 
while the rectangular nanoparticle showed a ^ conver- 
sion efficiency, as one might expect due to the absence of 
symmetry breaking. 

As in Ref.El] we find that an array of U-shaped gold 
nanoparticles can enhance the SHG efficiency by about 
two orders of magnitude with respect to other non- 
centrosymmetric nanoparticles, such as T-shaped parti- 
cles for example. Moreover, the polarization state of the 
generated field is rotated by 90° for a pumping electric 
field polarized along the x-direction. This polarization 
flip is not surprising, because symmetry breaking occurs 
along the y-direction. However, the nature of the con- 
version enhancement remains an unsettled point in the 
literature. In Ref.!^ the authors speculated that this ef- 
fect arose as a result of a strong nonlinear magnetic re- 
sponse associated with the Lorenz force. However, the 
Lorentz force, proportional to Pi x Bi, is but one of sev- 
eral influential contributions, Eq. ([3|. Within the context 
of our model, by simply switching off the Lorentz force 
contribution to the generated field we find that the con- 
version efficiency is only marginally affected. The effect 
is shown in Fig. [2^, where the second-harmonic conver- 
sion efficiency is calculated as a function of the pump- 
ing wavelength, Aff, with and without the Lorentz force 
term. Therefore, the magnetic response of the U-shaped 
nanoparticle does not seem to play an important role in 
the description of the SHG process, at least assuming the 
hydrodynamic model. 

The origin of the conversion enhancement may be ex- 
plained by studying the nonlinear electrical properties of 
the nanoparticle. If the Lorentz contribution is neglected, 
all the electric nonlinear contributions are basically pro- 
portional to the product of the electric field E and the 
induced charge distribution p = — en = V • P, as may 
easily be deduced from Eq. \2h\ . In Fig. two parallel 
situations are depicted: i) the rectangular nanoparticle 
and ii) the U-shaped nanoparticle. The arrows indicate 
the electric field at the surface, while the color map rep- 
resents the charge distribution. At the fundamental fre- 
quency both structures behave similarly. The external 
electric field of the electromagnetic radiation produces a 
gradient in the electrical potential that causes charges 
to migrate from one pole to the other. To a first ap- 
proximation, the resulting second-harmonic polarization 
is proportional to the product between the charge distri- 
bution and the electric field, that is: P(2a;) cx p(a;)E(a;). 
As shown in the sketch of Fig. [2]3, the nonlinear polar- 
ization will result in a second-harmonic electric field that 
cancels out for the centrosymmetric nanoparticle, result- 
ing in a near-zero second-harmonic efficiency. On the 
other hand, for the U-shaped nanoparticle the resulting 
harmonic fields oscillate in phase and the produced field 
can radiate to the far-field. These findings lead us to sug- 
gest that the characteristic emission of nonlinear optical 
SRRs is driven by the electric properties of the structure 
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FIG. 2. (Color on line) (a) Comparison of the second- 
harmonic conversion efficiency with and without Lorentz force 
contributions for the U-shaped nanoparticle of Fig. [l] (b) 
Sketch of the fundamental and second-harmonic currents flow- 
ing in the rectangular and U-shaped nanoparticles respec- 
tively. For the symmetric nanoparticle the nonlinear currents 
cancel out and no second-harmonic field is produced. Instead, 
for the U-shaped nanoparticle the currents oscillate in phase 
and radiate to the far-field. 



rather than its magnetic response, as originally conjec- 
tured. 

In conclusion, we have discussed the influence of the 
electron gas pressure on SHG and find a new way to 
express nonlinear surface contributions in terms of the 
polarization vector in bulk regions. Using this approach 
we have analyzed the SHG in several types of 3D gold 
nanoparticles and demonstrated that the origin of the 
previously experimentally reported enhanced SHG from 
SRR-based nano-structures is mostly electrically driven. 
This development should help simplify the investigation 
of the SHG process in full 3D metal structures, and en- 
able the investigation of arbitrary nano-structured plas- 
monic media, offering a viable design approach to inte- 
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grated, efficient nonlinear optical media. 
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